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AlWTHACT 


Ifii  vwMkwrttioM  with  thr  application  of  linear  chirp 
flllr#*  to  <ft*per*r«l  surface  wave  trains  recorded  at 

M«wata*»a,  a  nutnlN'r  of  eventa  from  the  SinoSoviet 
Ha%r  Wrti  r%atttined  in  order  to  determine  the  pre- 
•HilaMlili  of  thi-  time  of  the  maximum  filler  output.  Th(‘ 
^rocip  %r|o(itie*  thu*  oMained  were  compared  with  the 
tf^roimp  «rt(.MrStie*  deterodned  hy  an  analyal  using  ordinary 
r*ion  analv*ift. 

<lw  the  lM*la  of  he  scatter  in  group  velocities  found, 
**arriial  time  windows"  for  surface  waves  from  various 
ft  ti(Mi«  are  jiven 

l)i*p«^rsion  analysis  has  shown  that  the  assumption  of  a 
lifiear  frerfuem  y  variation  in  surface  waves  is  not  always  sa- 
ll•faltory-  Therefore,  we  have  allowed  for  non-linear  frc- 
*|Mew«  variation  in  the  chirp  filters,  and  some  applications 
c*f  •uch  filters  are  demonstrated. 
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|i  !♦  <Mr||  |i^»l  !*»'  1^  frl*IHr  r^l* 

^  r»<Ml  »*%r*  I*  #  |i>»'^vrf» 

f!)il  4l*l  f'it^l*«*M  l>w-|#rrf^  0^$ 

OsrJ*  li  I*  l«l♦w  I***!  ni'flrr  i*#*" 

**#%%  ^«lr  fto  rii^  I*  *  (#f$ri«<> 

I*  *h^  t0^h^l^  r  *4  iHr 

littiliitdv  l«  iM*  ♦iMc  rSfftlf«»fH  I*  l^r  drirt  Irwl 

•4  ih«*  •MfUfr 

In  (Hi*  i*orK  nr  ht«%r  I^c4ir4  Into  *(*fnr  at  lHr>  f»f#l4rm* 
conn-'Clril  wiiH  lh«*  •li'lrfilinn  i4  lH»*  n^vr*  THi* 

«kark  foti&Mk*  nionie  lin<*s  flr^lly  n  flifrM**inin  ul  (hr 
po**fHiUiy  of  prrrficline  (Hr  «rrl%°nl  limr*  of  *Mrf;irr 
uAvr  irAin*  from  »f»relfic  f*4M*  of  (hr  oorld:  nnH  *rrooHly 
(hr  rffrcl  of  Applying  dlffr rrnlly  conifirucird  chirp  fill<*r» 
for  different  reglun*  ua»  InveftllpAird. 
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II,  m  Till  i^Hi  im^i  AiMi*rTY  or  tiii.  okoui* 

VI  uoeri  V  oi  mi  ioirc  havi.moii  wavu 

Hi  r«ififir(tiiMi  iHr  «»|»|ilit!f|f<«»i  <«|  tiji«|f  Ht'il  filirr* 
l#»  *«irNrr  »>!»%r  <il  iHr  l^^ri^*  Af«rrl«irr 

({.AliA),  M««*«l4*i4.  <*  liMtttlurr  tt4  rvrti>i»  frotti 
III*"  iifiiliovlwl  irfilMi  l»?»vr  H<if*‘|i  rMitftiltirrf  In  utHrr  lO  «lr|rf  ■ 
ifilnr  ili*"  firriliiiiaMlili  il**'  ilwMf^  *4  iH*^  filirr 

Clr«iirl|-*  ll  ll«r  filirr  in*lfr4l  !♦<»•  •rvrrni  r*’« 

*M«llli»t  for  rv*tt'|»lr  frotti  Inirrfrrinf  rvrM*,  fl  Ir  f|r*irAUIr 
l»  Hr  <iiMr  lo  fr*irtci  III*"  r»|irrir«|  pr^l^  frotit  A  p^rllcuUr 
rvrM  lo  <1  •m«ll  rrifion  on  llir  littir  «*!l* 

Tllf«  filirr*  orr  #|r»||^lir||  |«  tmicH  *lf*t»r  r*r^  *urf<tcr 
ir«*itt»«  A**wttilni  tliirjir  grotip  d*"l<iiv  In  ihr  frr*)itrncv 
riintr  0.  02^  •  0,0%S  II#  (I)  THr  filirrlng  I*  prrfurtnrd 
In  Ihr  frr<|Mrnc|*  ilottoiln.  *nrf  Ihr  «rro  drUy  frrr|urncy  u*rd 
I*  0  0^  lit  Thr  filirr  Irngih  rrfrrrrd  lo  l*lrr  In  Ihl*  Irsi 
I*  Ihr  Irnglh  of  Ihr  limr  IttifniUr  rr»fH*n*r  of  Ihr  filirr  In  *rcond 
rtir4»t«rrd  from  Ihr  It  »rc  lo  Ihr  40  »rc  w^rr.  Thr  group 
vrioclly  found  hy  urSf  ihr  mtilmum  filirr  oulpui.  corrr* 
upond*  lo  Ihr  group  vriorlly  found  by  tn  nn^lytl  Idrnlifying 
Ihr  nrtivtl  limr  of  Ihr  20  »r(  fundtmrnl«l  modr  Htylrlgh 
w«vr»  on  Ihr  rtw  Ir^cr*  or  Ihr  brtm.  Tho  long  prrird 
brtm  Ift  formrd  by  uulng  n  phttr  vrioclly  of  1.  7  km  »rc 
«cro*»  Ihr  nrrny  for  nil  frr<iurncirii. 

In  Ihif  comm  Ihrrr  nrr  crrinln  qurfllon*  lo  br  nntwrrrd: 

I.  Civrn  iwo  or  morr  rpicrnlrnl  rrglont,  whni  in  ihr 

diffrrrncr  in  ihr  group  vrloclilr*  briurrn  Ihr  rrglons  nnd 
Ihr  army,  comparrd  lo  Ihr  ncailrr  in  group  vrlucilira 
wiihin  rach  rrgion? 


«  |«  |irrci*ii4M  wilh  wMrH  lh»'  llin**  of  iho 

fll^rr  Ci*M  I*  for  n  fMMirulAr  roKiott? 

i  ll*r  lo  I  nM  £,  I#  li  worllt  wltll««  lo 

world  «»•  •rf^n  ff«HM  iHi*  ^rray  in  Irrinii  of  Kroup 

%rtof 

4  Hflirrr  dor*  |lt««  filipf  coiop.’irpd  lo  lln*  flriti 
■♦rri^lt*#  iO  rr#  nnvfPt  o*  picked  liy  nn  nnnIyMi? 


IIA  t  \ 

\  folirfllon  t«f  rvrni*  frooi  Cpnir^l  Astl.i,  iltp  Kurile 
|»|*f»d*  fPttlon  i»nd  iHr  •oulhHp*irrn  liortler  count rIeM  of 
llsr  I'niofi  lr»*  ItTm  r»«»*lderrd.  I.*c«»pi  for  noitie 

P4fdo*|«m*  ff»*#ft  r4*|pr«  KoMlilt.  oi>l/  event#  for  which 
lltr  C09*l  4  C^eodellr  Survey  (USCGS)  epicenter* 
wrrr  IiMowm.  l♦o%=e  l*rrn  u*rd.  No  pov  fly  recorded  or  deep 
rvrinl  (I*  XMd  kitit  w**  included  in  the  collection.  The 
rvrni*  were  divided  Into  three  reition#  nccordinj;  to  the 
type  of  W  4>%rp4th*  to  I^SA  A*#utoin|:  cotttinentni  structure 
hene*iih  the  Artie  fkeon  the  p**h#  <ire: 


\  Crnlf.«l  A*i« 

H  Kurile  |*|«nd# 

C  Sonihw*  •tern  border 
«o<*n1fie*  of  U$SK 


Wrtveoftth 
ContinentAl 
l*rtrtly  ocennic 
North  Pticific  Oce.m 
l*Artly  ocennic 
North  Atinntic  Ocean 


A  Noinya  /.efolyn  and  a  weniern  Kazakh  explosion  have 
hern  included  in  reeion  A. 

'The  event*  are  li*ted  In  Tahle  I. 


PROCEDURE  AND  RESULTS. 


To  find  the  arrival  times  of  the  20  sec  waves  from  the 
raw  traces  or  the  beam,  the  standard  method  of  Ewing  and 
Press  has  been  used.  Arrival  times  of  successive  crests  are 
plotted  versus  crest  number,  and  the  resulting  points 
approximated  by  a  smooth  curve.  A  straight  line  with  a 
slope  corresponding  to  20  sec  period  is  then  adjusted  until 
it  is  tangent  to  the  smoothed  data.  The  point  of  tangency 
gives  the  arrival  time  of  the  20  sec  wave.  In  many  cases, 
once  the  period  of  the  surface  wave  train  lecorded  at  LASA 
has  dropped  to  20  sec,  long  portions  of  the  smoothed  curve 
tend  to  be  parallel  to  the  "20  sec  slope  line".  This  is  of 
course  most  pronounced  for  wave-paths  with  substantial 
amounts  of  oceanic  crustal  structure,  because  of  the  rapid 
dispersion  for  waves  of  about  20  sec  period  (Fig.  1). 

The  arrival  times  picked  by  the  analyst  represent  the 
beginning  of  these  waves.  Thus  one  can  expect  the  group 
velocities  found  by  an  analyst  to  be  slightly  larger  than 
the  group  velocities  calculated, using  the  time  of  the 
matched  filter  output.  USCG3  epicenters  and  origin 
times  have  been  used  in  calculating  the  group  velocities. 

The  average  and  the  extreme  values  of  group  velocities 
and  the  chirp  filter  lengths  for  the  different  regions  are 
given  below.  The  explosions  and  earthquakes  with  continen¬ 
tal  paths  have  been  separated  for  comparison. 
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CONTINENTAL  PATHS 


Group  velocities  km/sec 


Matched  Filter  Analyst 

Filter 

length (sec) 

High 

2.  96 

3.  04 

700 

Central  Asian 

Average 

2.  93 

3.  00 

560 

earthquakes 

Low 

2.  90 

2.  93 

350 

5  events 

High 

2.  99 

3.  13 

650 

Explosions 

Average 

2.  95 

3.  01 

505 

9  events 

Low 

2.  90 

2,  96 

350 

High 

2.  99 

3.  13 

700 

Earthquakes  and 

Average 

2,  94 

3.  00 

525 

explosions 

Low 

2.  90 

2.  96 

350 

14  events 

PARTLY 

OCEANIC  PATHS 

PACIFIC  OCEAN 

High 

3.  29 

3.  34 

700 

Kurile  Islands 

Average 

3.  15 

3.  17 

415 

earthquakes 

Low 

2.  97 

3.  01 

250 

16  events 
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NORTH  ATLANTIC  OCEAN 


High 

3.  11 

3.  16 

800 

• 

Southwestern  USSR 

Average 

3. 02 

3.  07 

535 

border  countries  earth 

Low 

2.  89 

3.01 

400 

quakes,  10  events 

COMMENTS 
Region  A. 

It  is  observed  that  the  scatter  in  group  velocities  be¬ 
tween  the  explosions  is  slightly  larger  than  for  the  earth¬ 
quakes.  However,  the  latter  group  consists  of  only  5  events, 
and  more  data  may  change  this.  At  a  typical  epicentral  dis¬ 
tance  of  85°,  the  matched  filter  peak  output  would  be  con¬ 
fined  to  a  time  window  of  100  sec  width,  for  the  group  veloci¬ 
ties  shown  above.  Thus  by  using  the  group  velocity  2.  94 
km/sec,  it  should  be  possible  to  predict  the  peak  output  within 
roughly  -  1  min  for  this  region.  The  results  support  the  idea 
of  a  largely  continental  crustal  structure  beneath  the  Arc'ic 
basin. 

The  (analyst)  velocity  3.  13  km/sec  for  one  of  the  Kazakh 
explosions  is  anomalously  high,  the  nearest  being  3.  05  km/sec. 

Region  B 

The  group  velocities  for  region  B  are  more  widely 
scattered.  This  is  the  case  for  the  hand-picked  group  ve¬ 
locities,  as  well  as  those  using  the  time  of  the  chirp  filter 
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output.  This  might  have  been  anticipated  because  of  the 
steepness  of  the  oceanic  dispersion  curve  (FiR*  !)• 

For  some  of  these  events  the  filler  outputs  show  multiple 
peaks  of  nearly  the  same  order  of  magnitude  as  the  main 
peak.  This  happens  to  be  the  case  for  the  events  with  the 
extreme  values  in  region  U.  Using  the  secondary  peaks 
for  these  events,  the  resulting  group  velocities  would  be 
close  to  the  average  (3.  I  *>  km/sec).  However,  the  extreme 
values  found  by  the  analyst  represent  th»‘  same  events,  so 
no  support  is  given  to  the  idea  of  restricting  the  "arrival 
time  window"  this  way.  Consequently,  one  has  to  expect 
this  scatter  in  group  velocities,  which  for  a  typical  diuiance 
of  65°  corresponds  to  4  min  on  the  arrival  tinte  axis 
(-  2|  min  using  the  average  velocity  given). 

Region  C. 

The  oceanic  parts  of  the  wavepalhs  are  relatively 
smaller  than  for  region  II,  and  large  parts  of  the  crust 
under  the  North  Atlantic  Ocean  are  less  oceanic  than  the 
crust  under  the  North  Pacific  Ocean. 

The  scatter  in  group  velocities  correspwods  to  nearly 
4  min  at  the  distances  involved  (85®).  The  anontalous  low 
velocity  2.89  km  sec  corresponds  to  the  event  with  thr 
maximum  filter  length.  (07.04.68  03:42:00  4,  Aegean  Sea). 
On  removing  this  event  the  width  of  the  arrival-time  region 
drops  to  about  140  sec  (  -  1^  min),  and  the  average  filter 
length  to  505  sec. 
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Ill  ON  TMK  APPLICATION  OF  EXPERIMENTALLY 


CLNEUATEI)  CHIRP  FILTERS 


In  lh»*  pr»*vlouH  Hcction  we  suggested  mapping  of  the 
wtirld  in  lermM  of  group  velocities. 

VVe  proceed  by  demonstrating  the  use  of  this  information 
in  (lelertion  methods.  In  this  section  we  will  discuss  the  use 
of  chirp  filters  designed  to  match  surface  waves.  The 
filters  are  generated  on  the  basis  of  the  experimentally 
found  frequency  variations  in  the  wave  trains  to  be  matched. 

Capon  and  Green  (2)  and  Capon  and  al.  (7)  used  chirp 
signals  with  linear  frequency-sweep  waveforms  as  matched 
filters  for  dispersed  wave  trains.  They  found  that  these 
waveforms  give  a  satisfactory  match  to  the  surface  waves 
and  demonstrated  considerable  enchancement  in  signal-to- 
noise  ratio.  However,  dispersion  analysis  of  wave  trains 
from  various  areas  his  shown  that  the  linear  frequency 
assumption  can  noi  be  retained  in  all  cases.  Accordingly, 
in  the  design  of  chirp  filters  we  have  allowed  for  non¬ 
linear  frequency  variations. 


METHOD 

Capon  and  Green  (2)  used  the  reference  waveform 


fl)  R(t)=«< 


f,  -f  , 

sin  (  2  T  (f  t  +  °  t^)  1  ,  0<  t  <  L 

O  c  Lj 


0,  otherwise. 
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where  are  the  initial  and  final  frequencies  of  the 

wave  train  whose  time  duration  is  assumed  to  be  L  seconds. 


To  allow  for  non-  linear  frequency  variation,  formula 

(1)  should  be  modified  to 
r  t 

sin  (Ztt  J  F  (t  )  dr  ),  0<t<L 

0 


(2)  R(t)  -  J 


^0,  otherwise, 

where,  for  the  special  case  of  a  linear  frequency  variation 
we  have 


F(r)  =  f„  + 


‘r‘o 


F  (  T  )  is  the  frequency  as  a  function  of  time  in  the  chirp 
signal.  This  frequency  variation  is  found  by  conventional 
dispersion  analysis  of  wave  trains  from  the  region  for 
which  the  filter  is  to  he  designed. 

However,  the  curves  thus  found,  are  encumbered  with 
errors  partly  due  to  inaccuracies  in  reading  the  peak  arri¬ 
val  times  and  partly  resulting  from  the  procedure  used  to 
obtain  the  frequency  functions.  Corrected  frequency  curves 
can  be  formed  by  adjusting  the  chirp  signals  on  the  basis  of 
a  visual  comparison  with  the  wave  trains  they  are  supposed 
to  match.  This  has  been  done  whenever  it  has  been  found 
necessary.  Furthermore,  because  of  disturbances  in  the 
wave  trains  to  be  analyzed,  like  noise  or  weak  recordings, 
the  chirp  signals  found  often  cover  only  part  of  the  frequency 
range  in  the  wave  trains.  In  such  cases  the  frequency  curves 
should  be  extended  in  some  way,  and  this  has  been  done  by 
some  trial  extrapolations  of  the  curves,  preceding  a  visual 
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comparison  of  the  resulting  filter  and  the  actually  recorded 
wave  trains. 

Frequency  curves  for  several  regions  have  been  obtained 
this  way,  and  some  of  these  curves  are  shown  in  Figs.  2 
and  3.  As  can  be  seen  from  Fig.  2  the  assumption  of  a 
linear  frequency  variation  would  give  a  poor  fit  to  wave 
trains  coming  from  some  regions.  Fig.  3  shows  the  curves 
used  in  this  paper.  The  dotted  part  of  the  Rayleigh  wave 
curve  for  West  Pakistan  is  the  extension  of  the  original 
curve.  As  an  example  the  Rayleigh  wave  chirp  filter  for 
West  Pakistan  is  shown  on  the  figure,  and  for  comparison 
both  the  original  and  the  extaided  version  are  shown. 

DATA  AND  RESULTS 

The  data  used  in  designing  and  testing  the  chirp  filters 
have  been  taken  from  Kongsberg  Standard  Station.  Since 
these  data  are  on  paper,  we  had  to  digitize  the  seismograms 
and  this  has  of  course  been  a  limiting  factor  on  the  data  base 
available,  and  accordingly,  on  the  testing  possible  to  perform. 

In  Table  2  the  events  used  in  this  study  are  listed.  Origin 
times  and  epicenter  positions  were  taken  from  the  preliminary 
lists  of  USCGS,  except  for  some  of  the  explosions  where  these 
data  were  taken  from  the  bulletins  of  the  International  Seismo- 
logical  Centre  ^ISC),  Edinburgh. 

4  different  epicentral  regions  have  been  selected,  of 
which  West  Pakistan  is  the  first  to  be  discussed.  From  this 
region  the  Kongsberg  recordings  cover  a  relatively  wide  magni¬ 
tude  scale;  the  events  of  13.  02.  66  and  08.  12.  66  are  both  given 
an  m^  =  5.  1  which  is  far  below  the  magnitude  of  events  nor¬ 
mally  detected  by  the  Kongsberg  long  period  system  at  this 
distance.  The  Rayleigh  wave  chirp  filter  for  this  region  is 
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based  on  the  event  of  24.  01. 66,  and  the  wave  trains  and 
filter  outputs  are  shown  in  Figs.  4-7.  The  amplitude  of 
the  surface  waves  from  the  07.  02.  66  event  is  compressed 
by  a  factor  of  about  3. 

All  outputs  are  given  with  an  output  frequency  of  0.  033 
Hz,  that  is,  the  peak  (or  through)  on  the  output  trace  is  lo¬ 
cated  at  the  arrival  of  the  30  sec  wave. 

In  the  ideal  case  the  output  of  the  filter  should  look 
like  an  autocorrelation  function,  and  we  can  see  that  this 
is  nearly  the  case  for  all  events  except  for  the  event  of 
08.  12.  66.  Despite  the  fact  that  this  event  is  given  the  same 
body  wave  magnitude  as  the  event  of  13.  02.  66,  the  surface 
waves  at  Kongsberg  are  much  weaker  and  are  in  fact  too  weak 
to  contain  very  much  energy  in  the  low  frequency  range  co¬ 
vered  by  the  filter.  To  detect  events  of  low  magnitude,  the 
chirp  filter  should  therefore  be  extended  towards  higher  fre¬ 
quencies.  The  extended  frequency  curve  and  the  resulting 
filter  are  shown  in  Fig.  3,  and  this  chirp  filter  matches  the 
surface  waves  from  all  three  events  fairly  well  (Figs.  8-10), 
although  the  high  frequency  content  in  the  chirp  signal  has 
modified  the  true  autocorrelation  form  of  the  output. 

The  other  epicenter  regions  selected  are  the  Nevada 
Test  Site  (NTS),  the  Rat  Islands  region  and  the  Nepal-India 
border  region.  The  frequency  curves  are  shown  in  Fig.  3, 
while  the  raw  traces  of  the  events  and  the  filter  outputs  are 
given  in  Figs.  11-18.  The  Nevada  events  are  all  underground 
nuclear  explosions,  while  for  the  Rat  Islands  region  there  are 
two  earthquakes  and  one  explosion.  A  Love  wave  chirp  filter 
has  been  used  for  the  Nepal-India  border  region,  Rayleigh 
wave  filters  for  the  other  regions. 
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In  each  of  these  cases  the  match  between  the  filters 
and  the  surface  waves  is  quite  good.  However,  the  events 
from  each  region  do  not  vary  very  much  with  regard  to 
magnitudes  and  positions  of  the  epicenters.  Therefore,  very 
little  information  can  be  drawn  from  these  examples  about 
the  consequence  of  using  the  same  chirp  filter  for  events 
which  differ  in  magnitude  or  epicenter  position.  As  to  the 
magnitude  variations,  however,  the  results  for  the  West 
Pakistan  region  as  expected  indicate  that  such  variations 
have  little  influence  on  the  output,  if  the  chirp  filter  is  de¬ 
signed  to  cover  a  broad  enough  frequency  range. 

The  azimuthal  range  of  each  region  possible  to  cover  by 
one  filter  is  of  course  dependent  on  the  crustal  structure  be¬ 
tween  the  epicenter  and  the  station,  and  must  be  tested  in 
each  case.  In  Fig.  19  the  results  of  applying  a  Love  wave 
chirp  filter  to  one  event  from  Northern  Italy  and  to  one  event 
from  Yugoslavia  are  shown.  The  filter  is  designed  for  North¬ 
ern  Italy.  The  two  events  have  about  the  same  epicentral 
distance,  and  the  result  shows  that  the  difference  in  azimuth 
does  not  seem  to  have  any  great  influence  on  the  wave  form 
in  this  case.  However,  in  cases  where  multipath  propa¬ 
gation  can  be  expected,  small  differences  in  azimuth  may 
change  the  waveforms  considerably  (8).  Among  the  wave 
trains  considered  in  this  paper,  the  wave  trains  from  the  Rat 
Islands  region  seem  to  consist  of  groups  of  waves  coming 
from  different  azimuths.  This  is  indicated  both  by  the  beating 
of  the. envelope  of  the  raw  traces,  and  by  the  form  of  the 
filtered  traces, which  shows  several  peaks  (Figs.  14-16). 

In  the  first  part  of  this  report  we  looked  into  the  pre¬ 
dictability  of  the  20  sec  Rayleigh  waves  from  different  epi¬ 
central  regions,  as  seen  from  LASA.  It  might  be  of  some 
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interest  to  calculate  group  velocities  for  different  paths  to 
Kongsberg  on  the  basin  of  the  chirp  filtered  events  where 
maximum  peak  (or  trough)  indicates  the  arrival  time  of 
the  30  sec  wave.  This  has  been  done  for  all  events  in 
Table  2  and  the  results  are  listed  to  the  right  of  the  table. 

The  results  show  very  little  scatter  in  the  group  ve¬ 
locities  within  each  region,  the  largest  variation  is  found 
for  West  Pakistan,  where  the  group  velocity  ranges  from 
3.  03  km/sec  to  3.  10  km/sec.  The  extremes  correspond 
to  the  events  of  lowest  and  highest  magnitude  respectively 
and  represent  a  difference  of  about  1  min  in  arrival  times. 

From  one  region  to  another  there  is,  however,  a  con¬ 
siderable  variation.  The  30  sec  Rayleigh  wave  group  ve¬ 
locity  from  West  Pakistan  is  about  3.05  km/sec  -  whereas 
the  corresponding  velocity  for  the  Rat  Islands  events  is 
about  3.  60  km/sec.  The  Love  wave  paths  investigated  in 
this  paper  are  all  continental  and  the  group  velocities  for 
the  30  sec  Love  wave  are  close  to  3.  5  km/sec. 
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IV  COHCUUSIOMI 


Chirp  fitirrinc  «»f*{»Ur<)  to  nrr^y  diilii  or 

•lw#lr>  •rH»or  «>»ilfMt*  In  both  €<*•»♦*  ih»*  prt^ituiord 
«*rrlv*l  t4  thr  •orl^rt^  o4Vi«*  U  of  jeri*nt  Itopor- 
l*nc**  tartly  birf<iiM*r'  or  tho*  <i»n  rKfludr  Inirrfrrinft 
r%rtft*«  And  partly  b(CAti*r  of  rAvIne  compulinfs  lintr. 
if*fn4f  ir<»fn»  frtordrd  At  t^SA  from  diffrrrni 
rfOrmt#'*!  fr^ien*  it  h»»  l*r«'n  foaml  ihAt  ihr  Arrival 
tiftir  of  thr  ih  •rr  oovr#  routd  hr  prrdicird 

within  roiifhly  •  I  ffiln  In  thr  hr*t  fA»r  And  •  2  min 
In  thr  wor*t  c#»r  thrrr  notnhrrA  Arr  of  foornr  In- 
flmnc 'd  hf  thr  *|#r  of  aM  dlrMncr  to  thr  rpIrrntrAl 
irfim*  And  #l-»o  by  thr  «OHtplr*lty  of  thr  #truciurr 
loiwrrfi  thr  •««<f»rr  And  thr  rrcrirrr. 

Ihr  fr*nlt«  IndtCAtr  thAt  thr  otr  of  non  llnrAr 
ff\r«ii'wnc%**A*rp  wA%^r  form*  a*  toAtchrd  flltrrA  con 
br  *9i  rffrrtNr  tond  In  drircilon  of  rorlACr  oAvrn. 

Il  ml#hl  W  obprrtfd  thrt  onr  rijOAtly  writ  cootd  u»r  thr  ACtu- 
Ally  frrofdrd  WA%r^  tf  AlO*  Aa  tOAtchrd  flltrrA.  thtt*  aI*o 

thr  AfwplIfHtidr  AAf  lAtlon*  lof  Inch'd.  Hoorvrr.  Cnpon 
And  aI  (Tf  rhoArd  th*t  AfofdllAidr  %‘AflAllenA  hAvr  vrry 
iHlIr  rfirfi  on  thr  io<nol*r  fAtlo. 

0h^<tlon*  rAl»"«%Ani  to  owr  tnAtfhrd  flttrrin|i  proc'^durr 
hAir>  hrrn  pd  forward  frtrHly  by  CApon  (1)  In  A  p^prr 
€  AntrfAisnt  tV  pAyl*l#h.<A>%o  rnnhlpAih  propAjtAilon  At 

f  h»  WA%A  trAl*!i*  ^  AnA|y#rd  waa  Ahooti  to  rontlAt 
of  ♦rAA  fAl  ###Ap*  O!  wA%r*  AfflAint  frtoti  thffrrrnt  a#J- 
fAjAtMii  dffr«li«A»*  thl*  fAHfiltip*th  pfopAfAtlon  hr  rAplAlnrd 
to  l>r  dr*  to  th*  fh**r'»vrlOf  Wy  fontf  AAl  brtwrro  ocrAnIc 


and  continental  regions  which  led  to  refractions  and  re¬ 
flections  at  the  continental  margins.  Also,  he  found  that 
if  the  multipath  propagation  for  a  given  epicenter  was 
known,  it  would  still  be  difficult  to  predict  this  effect  for 
an  epicenter  whose  location  was  significantly  different 
from  that  of  the  given  epicenter. 

These  results  make  the  use  of  master  events  as  matched 
filters  for  Rayleigh  waves  rather  doubtful,  as  emphasiz.ed  by 
Capon,  and  he  concluded  that  the  use  of  chirp  wave  forms 
for  matched  filtering,  as  recommended  by  Capon  and  al.  (7), 
may  represent  the  most  reasonable  approach  to  the  problem 
of  matched  filtering  of  surface  waves. 

The  objections  against  the  use  of  master  events  as 
matched  filters  apply  to  our  non-linear  chirp  filters  also, 
since  we  are  using  master  events  to  find  the  frequency 
curves.  However,  if  our  chirp  filters  are  to  be  used  on 
recordings  from  Kongsberg  or  the  Norwegian  Seismic  Array 
(NORSAR),  the  sites  of  these  stations  relative  to  the  epi¬ 
center  regions  of  main  interest  remove  at  least  parts  of  the 
objections.  This  is  due  to  the  fact  that  if  one  is  looking  for 
events  in  Eurasia,  the  paths  to  Kongsberg,  or  NORSAR, 
are  purely  continental.  Thus  no  severe  geological  ob¬ 
stacles  causing  multipath  propagation  should  be  assumed, 
and  this  is  crudely  corroborated  by  the  shapes  of  the  wave 
trains  from  these  regions  (Fig.  5)  in  contrast  to  wave 
trains  for  example  from  the  Aleutian  arc  (Fig.  14-16). 

Also,  in  designing  the  chirp  filters,  only  the  first 
Mucll-dcfincd  wave  packet  is  used  in  the  disperison  ana¬ 
lysis,  and  thereafter  the  filters  are  extended,  if  neces¬ 
sary,  by  extrapolation.  By  this  procedure  we  think  that 
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most  of  the  disadvantage  by  using  master  events  as  matched 
filters  is  avoided,  and  at  the  same  time  the  filters  have  a 
frequency  variation  with  a  better  fit  to  the  actual  wave  trains 
than  linear  frequency- sweep  wave  forms.  Furthermore  if 
we  want  tO  assemble  matched  filters  for  a  large  number  of 
regions,  actual  recordings  require  considerable  more  me¬ 
mory  space  than  our  frequency  curves,  which  can  be  re¬ 
constructed  from  a  few  points  (10-20)  on  each  curve. 
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Fig.  1.  Summary  of  observed  Rayleigh  wave  dispersion 
(After  J.  Oliver  1962  (6)  ), 


Tim*  (min.) 


Fig.  2,  Curves  showing  the  frequency  variation  found  in 

Love  (L)  and  Rayleigh  (R)  wave  trains  from  sc- 

0 

lected  epiccntral  regions. 
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Frequency  (Hz) 


Fig.  3.  Frequency  curves  used  In  constructing  the  chirp 
filters  applied  in  this  study.  The  original  and  the 
extended  filters  for  West  Pakistan  are  also  sho^^'U. 
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NONCSBCNC.  2N.OI.66.  LPI  WEST  PAKISTAN  N8«S.9 


TINCININ) 

Fig.  A.  Wmt  Piikistan  event.  Raw  trace  (bottom)  and  chirp 
niter  output  (top). 
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cnossconncLRTiON.  output  FncoucNcr  o.qsshz 


26 


flHill«iii«|  -triNi  •'• 

<fMc^  -MiM 


I  1C  ti:>  Wr*!  r%’t«rvt  iritcr  (bdllotn)  Anti  vh- 

ritirii  filtrt  (»>iit|»4t 


fl#  II  f  If *«>■  (♦M 

tin*  f  f!«^ 


29 


1  2 


8 


10  11  12  18  m  15  16  17  18  19  20 

TIHE  fMIN) 


KONGSBERC.  15.02.65.  LPZ  RRT  ISLANDS  HB>5.8 
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KBNCSBCRG.  02.10.69.  LPZ  riAT  ISLANDS  (CXPLQSIDNI  Nt>t.S 


TIACrAINI 


Fig.  16.  Rat  Islands  explosion  (Milrow).  Raw  trace  (Itultoitt) 
and  chirp  filter  output  (top). 


Fig.  Is.  Nrp^l'lftdWi  bordr  rrgtcm  lr-»cr 

(txictom)  404  chirp  ^if*rr  oafpitt  (t«ip) 
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CROSSCORnCLflTlON,  OUTRUT  FRCQUENCT  0.033  HZ 
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kROSSCORRELATlON.  OUTPUT  PREQUENCY  0.U33HZ 
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RONCSBERC.  30.  P.6. ,  LPEil  NORTHcRN  ITALY  46-  5.3 


1  I  5  C  S  b  ^  S  5  io  il  12  13  14 

TIME  IHIN) 

I  9.  An  example  of  applying  the  same  filter  to  two  events 
The  filter  has  been  constructed  on  the  basis  of  the 
Northern  Italy  event.  The  epicentral  distances  are 
approximately  the  same  while  the  azimuths  differ 
by  1 5  degrees. 


33 


